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ABSTRACT 

Knowledge of the microwave properties at cryogenic temperatures of components 
fabricated using High-Temperature-Superconductors (HTS) is useful in the design of HTS- 
based microwave circuits. Therefore, fast and reliable characterization techniques have been 
developed to study the aforementioned properties. In this paper, we discuss computer analysis 
techniques employed in the cryogenic characterization of HTS-based resonators. The revised 
data analysis process requires minimal user input, and organizes the data in a form that is 
easily accessible by the user for further examination. These programs retrieve data generated 
during the cryogenic characterization at microwave frequencies of HTS based resonators and 
use it to calculate parameters such as the loaded and unloaded quality factors (Q and Q 0 , 
respectively), the resonant frequency (f 0 ), and the coupling coefficient (k), which are 
important quantities in the evaluation of HTS resonators. While the data are also stored for 
further use, the programs allow the user to obtain a graphical representation of any of the 
measured parameters as a function of temperature soon after the completion of the cryogenic 
measurement cycle. Although these programs were developed to study planar HTS-based 
resonators operating in the reflection mode, they could also be used in the cryogenic 
characterization of two ports (i.e., reflection/transmission) resonators. 


INTRODUCTION 

The discovery of High-Temperature-Superconductors (HTS) in 1986 has resulted in 
passive planar microwave devices such as resonators and filters with higher quality factors 
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Figure I. YBCO microstrip linear, ring, and conductor-backed cophmar waveguide resonators 
deposited on lxl cm, 508 pm (20 mil) thick, dielectric substrates 


and lower insertion loss than those fabricated using normal conductors *" 3 . Significant 
attention has been given to the processing and characterization of HTS-based planar 
microwave components as efforts continue towards the inclusion of HTS-based circuits in 
working systems for space- and ground-based communications. Therefore, fast and accurate 
characterization and more advanced data analysis capabilities become extremely important in 
this endeavor. This paper deals with the data analysis programs used in our laboratory in the 
study of HTS-based microwave resonators. Extensive modifications of preexisting programs 
(written in-house for data analysis purposes) were performed in order to expedite the 
characterization of these resonators without sacrificing accuracy. The structure of these 
programs will be described, and data taken using YBa 2 Cu 3 () 7 _ A (YBCO) resonators on 
lanthanum aluminate (LaAl() 3 ) and R-plane sapphire (R-Al ? 0 3 ) substrates will be presented. 



Closed cycle 
helium gai 
refrigerator 

Figure 2, bxpcrimental setup used in cryogenic testing of HTS-based resonators. 
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Figure 3. Flowchart for data analysis program ’QANALYSIS'. 


EXPERIMENTAL 

Microstrip linear and ring resonators, as well as conductor-backed coplanar waveguide 
(CBCPW) resonators were fabricated on YBCO thin films deposited by laser ablation on 
LaAl() 3 substrates and on R-AI 2 0 3 substrates (with Ce() 2 buffer layers), respectively, using 
standard photolithography and chemical etching techniques. A schematic of these resonators 
is shown in Figure 1 . These resonators were tested by measuring the reflection scattering 
parameter (S,,)_ using the experimental configuration shown in Fig 2 4 . Briefly, the setup 
consists of an HP-8510C network analyzer, an HP 9000-216 computer, and a DRC-91C Lake 
Shore cryotronics temperature controller all coupled through the use of an IEEE-488 data bus, 
along with a closed-cycle-helium-gas refrigerator containing the device under test (DUT), 
which allows for cry ogenic measurements in the temperature range from 300K. to 20K. The 
data measured through the use of this system were then used by the data analysis programs to 
determine the loaded and unloaded quality factors (Q and Q„, respectively), the resonant 
frequency (f). and the coupling coefficient (k) of the resonator. 

Figure 3 shows a flowchart for the data analysis program we call 'Q ANALYSIS', after 
modification and automation. When automating this program, the first step was to have the 
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computer choose the files to be analyzed from the disk. In this way, the user would not have 
to be familiar with the format used in naming the files in order to choose those which are 
important in analyzing the data. Data files are named in the following way. The first four 
characters are always ’DATA'. The next 2 or 3 characters are numerals which show the 
temperature at which the data was taken. The next character is either a 'W T or a 'C\ showing 
whether the data were taken during the warming or cooling cycle. The final two characters 
are either ’1 V or '21\ indicating whether the data are S n (reflection) or S 21 (transmission). 
Thus, f Q ANALYSIS’ can be used for the analysis of one port (refection mode) as well as two 
ports (reflection/transmission modes) resonators. Since in this paper we are dealing with one 
port HTS-based resonators, ’Q ANALYSIS' catalogs the disk looking for files w'ith the first 
four characters 'DATA' and the last two 'll'. It counts all of the files which fit that 
requirement, and sets up a counting loop which will analyze only those particular files. 
When all of the appropriate files have been analyzed, the program ends. The program may 
terminate execution before all of these files have been analyzed, if certain conditions exist 
in the data, as will be described later in this paper. 

After selecting the data files corresponding to the particular temperatures of interest, 
the data for the first file in the array is displayed as a Smith Chart plot on the computer’s 
CRT screen. The Smith Chart’’ facilitates the graphical solution of transmission line 
problems and consists of a plot of the normalized impedance, given by 

Z/Z 0 = r + jx (1) 

w ; ith the magnitude and phase of a generalized complex reflection coefficient in a unity 
circle. The chart is useful for the analysis of lossless and lossy lines. It is used in this case 
to determine the coupling status of the resonator (i.e., overcoupled (k>1), undercoupled 
(k<I ), or critically coupled(K=l)) where k, the coupling coefficient, is defined as 

k = ( i±|r r |) / (i + |r r |> (2) 

with the sign depending on the coupling status of the circuit, f, is the reflection coefficient 
of the circuit at resonance r \ For k=I, T r =() indicating that the resonator is perfectly matched 
to the source. ’QANALYSIS’ calculates the radius of the circle corresponding to the 
experimental data, then compares that value with that of the circle on the Smith Chart 
corresponding to critical coupling, which has a radius of 0.5. If the radius of the plot of 
the experimental data is less than 0.5, the circuit is undercoupled, and if it is over 0.5, the 
circuit is overcoupled. Next, the program plots the magnitude of the data on a logarithmic 
scale as a function of frequency. Depending on the quality of the resonator, this plot 
should show a minimum at or very near the resonant frequency for which the resonator was 
designed. QANALYSIS' then goes on to determine f 0 . This step is implemented such that 
the correct resonance is selected; i.e., after finding the frequency at which a minimum point 
occurs, the program searches the entire frequency band under study for any other 
resonance. This accounts for cases in which a split resonance may occur, as for example, 
when testing microstrip ring resonators, where secondary or spurious resonances can appear 
due to variations in thickness, width, etc., of the line 7 . These secondary resonances could 
be larger than the primary resonance when the circuit is overcoupled. Note that with minor 
modifications, this searching capability could be applied in the analysis of multipole 
resonant circuits. The program also verifies that resonance still exits in the data; if the 
minimum point is not more than 3 dB below the value of the data where it is tlat (i.e., at 
a frequency far from f 0 ), then the program indicates that a resonance no longer exists and 
the analysis is completed. 

After the determination of f 0 is complete, it is then used to calculate the loaded quality 
factor of the resonator, since 
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Q ANALYSIS FOR FILE:DATA77Ul! 76. 7 KELVIN 
RESCNANT FREQUENCY IS; 7.9G629G19G7S 6H* 

UNLOADED 'O' IS: 14G0.35SG0716 

LOADED 'Q' IS: 449.5IS607B37 

COUPLING COEFFICIENT IS: 2.29368265B96 

COUPLING LOSS IS; .0997639596878 

FIT; HIGH 

HALF POUER SPAN IS: .W5381338B3643 . 01772J9S66395 

Figure 4? Sample printout of data from ’(^ANALYSiS'. 


Q = (f„ / At) (3) 

where Af is the half power frequency span. 

The final step in analyzing each data file is to fit a curve to the resonance. This 
assumes that the graph of the resonance resembles a parabola, and uses a second degree 
polynomial for the curve fitting. Depending on the shape of the curve, e.g. sharp, 
moderate, or shallow, the program will classify the fit as either High Q, Mid Q, or Low Q 
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Figure 5. Flowchart for plotting program 'DATAPLOT. 



respectively, which will be used to aid in the determination of the unloaded quality factor. 
In calculating Q„, the coupling effects must always be considered, regardless of the coupling 
state of the resonator. The coupling loss (a), which modifies the location of the half 
power frequency points, can be found by 

cj =(1-1^1) /(l + |rj) (4) 

where is the reflection coefficient of the circuit far from f t) In this program, Q 0 is 
calculated in a similar manner to Q, except that for the calculation of Q o , the half power 
frequency span is that obtained after correcting for coupling losses in the circuit. The 
program also calculates several other parameters, then prints a list which includes the file's 
name, the temperature, Q, Q u , f () , k, 0, the data fit type, and the half power span (3 dB 
points) used for determining Q 0 and Q, respectively, as seen in Figure 4. These parameters 
can be printed to the CRT or an external printer for further use. If there are more data files 
remaining, the process is repeated; if not, arrays of Q, Q t) , f, , and k versus temperature are 
stored to a floppy disk. The user is then given the option of plotting any of the stored 
parameters versus temperature or terminating the program. If the user wants to plot any 
of the stored quantities, then the program DATAPLOT' can be called to do so. 

Figure 5 is a flowchart of 'DATAPLOT. The user enters which of the parameters 
stored by ' QANALYSIS' is to be plotted (all are plotted as a function of temperature), the 
graph's title and legend, and whether the data is to be plotted to the CRT or to a plotter. 
Once the program has verified with the user that the information entered is correct, it loads 
the data from the chosen file and reads it. It determines the maxima and minima of the 
data set. then chooses a graphing range which will center the data on the graph, as well as 
making the graph easy to read. Finally, it plots the data and labels the graph with a title, 
legend, and any other user-specified comments. 

Figure 6 shows an actual plot from DATAPLOT 1 of Q 0 for both a linear and ring 
resonator of YBCO on LaA10 v Figure 7 shows plots of experimental data for a YBCO 
CBCPW resonator on R-ALO,. These plots give an idea of the format and scope of the 
plots produced by DATAPLOT 1 . The reproducibility of these data was determined by 
performing multiple measurement cycles. The reproducibility of Q tt was found to be within 
two percent at 77 K, and within 1 .5 percent at 20 K. The reproducibility from run-to-run for 
f 0 was below 0.02 percent at the aforementioned temperatures. 
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Figure 6. Sample data plot for the unloaded Q corresponding to YBCO linear and ring resonators. 
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Rcfoa&nt Frequency ; 6.9 GHz Resonant Frequency : 6,9 GHz 

Figure 7. Plots of experimental data for a YBCO conductor backed coplanar waveguide resonator 
(CBCPW). 


CONCLUSIONS 

In the course of this work, automated data analysis programs which increase efficiency 
during the characterization of superconducting microwave resonators were developed. These 
programs significantly reduce the time commitment of the user in the analysis process, and 
provide quick access to the results of the cryogenic characterization of HTS-based resonators. 
Although these programs were originally developed for the study of one port (i.e., reflection 
mode) HTS-based resonators, they could also be used for the cryogenic characterization of 
two port (i.e.. retlection/transmission) resonators. 
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